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Abstract 

We consider a SU(2)®U(1) gauge theory with many scalar doublets, but without 
fermions. We devise a systematic method of constructing quantities which are invariant 
j — , under changes of basis of the scalar fields, but acquire a minus sign under a general CP 

CN \ transformation. Those quantities provide a basis-independent way of checking for the 

■ existence of CP violation. 

Oh! 

p; 1 Introduction 

<D ' 

In 1985, Jarlskog M and, almost simultaneously, Dunietz, Greenberg, and Wu discovered 
that CP violation in the quark sector of the three-generation standard model originates in a 
^ ■ single basic CP-violating quantity, 

& : 

J = det[M u MlM d Ml] 

I 2 2\/ 2 2\/ 2 2\f 2 2\/ 2 2\/ 2 2\ 

oc [m t - m c )(m 4 - mj(m c - m u )(m fe - mj(m fe - m d ){m s - m d ) 

xIm(V ud V cs V* s V* d ) , (1) 

where M u and Md are the mass matrices of the up-type and down-type quarks, respectively, 
and V is the Cabibbo-Kobayashi-Maskawa matrix. Shortly afterwards, Bernabeu, Branco, 
and Gronau (BBG) found J by searching for quantities which, though invariant under a 
basis transformation of the fermion fields, acquire a minus sign under a CP transformation. 
Their method of deriving basic CP-odd quantities in the fermion-gauge sector of the standard 
model works for any number of fermion generations, though unfortunately for more than three 
generations a minimal set of independent CP-odd quantities is difficult to find ||. 

The methods of BBG have later been used in the construction of basic CP-odd quantities in 
models with vector-like quarks f|, with left-right symmetry 0, with neutrinos with Majorana 
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masses 0, and in supersymmetric models ||. With different aims, general CP transformations 
have also been extensively studied by Grimus and collaborators ||. 

Recently, Mendez and Pomarol |10| studied CP violation in the scalar sector of a SU(2)(g>U(l) 
model with many scalar multiplets. They only considered the possible clash between the CP 
transformation properties of the scalar gauge interactions and of the scalar mass matrices. In 
particular, for the simple case of the two-Higgs-doublet model, they phenomenologically derived 
the existence of a basic CP- violating quantity Ji, analogous to the J of Eq. [1]. 

In this work we consider the derivation of J\ as a basic CP-odd quantity under a general CP 
transformation, along the same line of reasoning by which BBG found J from basic principles. 
We further consider the CP violation originating not only in the mass matrices of the scalars, but 
also in the cubic and quartic scalar interactions, and do so for any number of scalar doublets. For 
simplicity, we only study a SU(2)®U(1) model with many scalar doublets, and do not consider 
the possibilities of either an extended gauge group, or of scalars in other representations of the 
gauge group. We also do not consider the Yukawa interactions of the scalars with the fermions, 
though in general these will constitute further sources of CP violation. This we do in order to 
keep our study within the limits of reasonable generality. However, our methods, supplemented 
by the ones of BBG, can, in principle, be applied to CP violation in specific models with both 
scalars and fermions. 

In section 2 we explain our method for the neutral sector of the two-Higgs-doublet model. 
In section 3 we extend the method to the multi-Higgs-doublet case. In section 4 we extend it 
to the charged sector. Section 5 summarizes our results. 



2 Two-doublet case 

We consider a SU(2)®U(1) gauge model with two scalar doublets Hi and H 2 . We do not impose 
any symmetry beyond the gauge symmetry. The most general Higgs potential consistent with 
renormalizability is 

V = ^HlHt+^H^ + i^Hl^ + h.c.) 

+\ 1 {h\h 1 ) 2 + \ 2 {h\h 2 ) 2 + x 3 (hIh 1 )(hIh 2 ) + x a {h\h 2 ){hIhi) 

+ [\ 5 (hIh 2 ) 2 + x & {h\Hi){h\h 2 ) + x 7 {hIh 2 ){h\h 2 ) + h.c] , (2) 

in which all the coupling constants, except /i 3 , A 5 , A 6 , and A7, are real by hermiticity. Without 
loss of generality, we assume that the two scalar doublets are in the basis in which only one of 
them, Hi, has a non-zero vacuum expectation value (VEV), v, which is real. Thus, 

Hl = {v + (H° G +iG°)/V2) ' (3) 
( (R + iI)/V2 ) ' (4) 



H 
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where H°, R and / are real neutral fields, and G + and G° are the Goldstone bosons, which, in 
the unitary gauge, become the longitudinal components of the W + and of the Z°. The stability 
conditions of the vacuum read 

fXi = -2Xiv\ (5) 
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fi 3 = -X 6 v 2 . (6) 

We use these conditions to eliminate Hi and fi 3 as independent variables from V . The fact that 
Eqs. [5] and |6] appear to overdetermine v is connected to the fact that we have chosen the basis 
in which H2 has vanishing VEV. This does not represent any loss of generality. 

Let us consider the potential V as a function of the fields £1 = R and £ 2 = J. V can be 
written as a sum in which some terms are of one of the following formsQ: 

1. aS, in which 5* denotes a general CP even monomial of H°, (G°G°), (G~G + ) and (H~H + ) 
(or a constant), and a denotes the respective coefficient, which is a polynomial of the 
coupling constants in V; 

2. bi^iS, terms which are linear in either R or J; the b{ are two coupling constants^; 

3. Cij£i£jS, terms quadratic in R and /, where the couplings c^- can be viewed as a symmetric 
2x2 matrix; 

4- dijk£i£j£kS \ where the couplings are symmetric under the interchange of any pair of 
indices; 

5. eijkiiiijthiiS ', with e ijk i totally symmetric. 

To be exact, the coefficients a through e^i should carry a superscript S specifying which 
monomial they refer to [for example, S = 1, or S = H°, or S = (H ) 2 , or S = H°H~ H + }. We 
will suppress that superscript in order not to clutter the notation, and will instead introduce 
primes in some specific examples. 

There is in V, for instance, a cubic interaction 

V2vH-H + [\ 3 H° + (Re\ 7 )R - (ImA 7 )/] , (7) 

which can be considered as the sum of a term of the form 1, with S = H°H~H + and a' = \/2v\3, 
and of a term of the form 2, with S = H~H + , and b[ = \/2vKe\j, b' 2 = — v^ImAy. 
The quadratic terms 

1 ( H °\ 

- ( H° R I ) M \R j , (8) 

where M is a 3 x 3 real symmetric matrix, are of particular interest, because they are the mass 
terms of the neutral scalars. By diagonalizing M, 

TMT T = diag(m?, m 2 2 , mj) , (9) 

where T is a matrix of SO (3), one obtains the squared masses m 2 of the three physical scalars 
X: 



fx, ^ 




( H 


x 2 


-T 


R 


\Xs J 







(10) 



1 Other terms, including the CP-odd field G°, such as G°-ff°£i, will be discussed later. Terms including the 
charged scalars, such as G~H + £ i i, will be considered in section 4, when we discuss the treatment of quantities 
involving the charged fields. 

2 We use throughout this paper the summation convention. 
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The mass terms in Eq. ||]can be viewed as the sum of a term of the form 1, (l/2)Mn(H ) 2 , a term 
of the form 2, H°(M 12 R + M 13 J), and a term of the form 3, (1/2) (M 22 R 2 + 2M 23 RI + M 33 J 2 ). 



The specific values of the coefficients are 

M n = 4X lV 2 , (11) 

M 22 = /ji 2 + (A 3 + A 4 + 2Re\ 5 )v 2 , (12) 

M 33 = fi 2 + (A 3 + A 4 - 2ReA 5 )t; 2 , (13) 

M 23 = -2 W 2 lmA 5 , (14) 

Mi 2 = 2w 2 ReA 6 , (15) 

M13 = -2v 2 ImA 6 . (16) 



Now consider a basis transformation. We do not want to leave the useful basis in which 
only Hi has a VEV, and that VEV is real; therefore, the most general basis transformation 
possible just rotates H 2 by means of a U(l) phase. This is equivalent to rotating R and / by 
means of a SO (2) matrix. Obviously, under such a basis transformation, the coefficient a of 
a term of the form 1 is invariant, while the two coefficients h\ and b 2 of a term of the form 2 
transform as a vector of SO (2). Similarly, the coefficients of a term of the form 3 transform 
as a second-rank tensor, the coefficients dijk of a term of the form 4 constitute a third-rank 
tensor, and so on. But the phase of H 2 is arbitrary and meaningless. Therefore, only quantities 
invariant under this basis transformation are meaningful and can appear in the final result for 
any physical quantity. For instance, a, or {pi) 2 + (b 2 ) 2 , or Cnc 22 — (ci 2 ) 2 , are invariant. 

Now consider a CP transformation. In the simplest CP transformation, besides the change 
of sign of the space coordinates, both Hi and H 2 transform to their hermitian conjugates^. This 
sets the CP properties for the Hi component fields: H° is CP-even, G° is CP-odd, G~ and G + 
are CP conjugates. For the H 2 components it dictates the transformations: H <-> H + , R —>■ R 
and / — > —I. However, because the phase of H 2 is meaningless, a general CP transformation 
must also include a phase rotation of H 2 . For instance, H <-> —H + , R — > —R and / — > I is 
also a CP transformation. In a general CP transformation, R and / get rotated by an 0(2) 
matrix with determinant —1. A CP transformation is analogous to a basis transformation, 
the only difference being that the 0(2) matrix which transforms the coefficients bi, Cij, and 
so on, has determinant —1 for a CP transformation, while it has determinant +1 for a basis 
transformation. 

Our purpose is to find quantities which are invariant under a basis transformation but 
change sign under a CP transformation. We need for that purpose the antisymmetric tensor 
eij, with ei2 = 1. This quantity behaves as a second-rank tensor under a basis transformation, 
but changes sign under a CP transformation: 

O ir O js e rs = (detO)eij . (17) 

It is easy to construct CP-odd invariants using e. For instance, bibje^, biC^b^^ e^jkCkieu, 
and so on, are quantities which, if they are non-zero, imply CP violation, while being invariant 
under a basis transformation. 

3 The hermitian conjugation arises in order to obtain the correct CP transformation of the gauge interactions 
of the scalars. In a model without gauge interactions, the hermitian conjugation would not be necessary, and CP 
conservation would be automatic. Indeed, CP loses its meaning if there are no gauge interactions to distinguish 
the particles from the anti-particles. 
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Let us apply this general method to the mass matrix M. As we pointed out before, that 
mass matrix can be divided in an invariant a, a vector foj, and a second-rank tensor Cy. The 
simplest CP-odd invariant which can be formed from this is 

Ji = b iCij b k e jk = b iCll b 2 - bicah = M 12 M 13 (M 22 - M 33 ) + M 23 [(M 13 ) 2 - (M 12 ) 2 ] . (18) 

This is the quantity introduced by Mendez and Pomarol ||10|| , who however did not explain its 
origin. J\ can be expressed as a function of the masses of the physical scalars and of the matrix 
T by using Eq. |. One obtains 

J x = mim^TiiTji(T & T j3 - T & T j2 ) (19) 
= (m\ - m\)(m\ - m\)(m\ - ml)T n T 21 T 31 . (20) 

Notice that, T being a matrix of SO(3), T i2 Tj 3 — T i3 Tj 2 = T fcl with (ijk) cyclic. 

The similarity of Eq. [2(] with Eq. [I] is evident. It is easy to show, using the identity 

8ik$ji — e ij e ki + Su5jk , (21) 

(Sij is the Kronecker symbol) that any other CP-odd invariant constructed from the mass matrix 
alone is proportional to J\. Thus, J\ occupies in the two-doublet model a place similar to the 
one of the quantity J of Eq. [I] in the three-generation standard model. If the masses of any two 
of the three physical scalars are equal, or if any of the matrix elements of the first column of T 
vanishes, there is no CP violation in the mixing of the neutral scalars. If any two of the three 
non-diagonal entries of the mass matrix M 12 , M 13 , and M 23 vanishes, there is no CP violation 
in the mixing. 

Having found that J\ ^ implies CP violation, the question arises whether J\ = implies 
CP conservation. In order to answer this, it is convenient to express J% as a function of the 
couplings in the potential in Eq. 0. From Eq. [18] and Eqs. [TT] to [16| one finds 



Jj = -8w 6 Im(A*A 2 ) . (22) 

The quantity in the right-hand side of this equation is invariant under a rephasing of H 2 and, 
being an imaginary part, it certainly implies CP violation. But one recognizes that there 
are other quantities with the same properties: the imaginary parts of (A5A7) and of (A^). 
Therefore, we guess that, even if J\ = implies CP conservation in the mixing of the neutral 
scalars, there will be other CP-violating quantities, appearing in the cubic and quartic scalar 
interactions, which are independent of J\. This point has not been emphasized before. 

Let us give a simple example. Consider the cubic interactions of the neutral scalars with 
the charged scalar given in Eq. [7[ As pointed out after that equation, those interactions consist 
of an invariant and a vector under a rotation of R and /. We may join the vector with the 
second-rank tensor from the mass matrix to form a CP-odd invariant: 

J 2 = V iCij V k e jk = b[b' 2 (M 22 - M 33 ) + M 23 [(6 2 ) 2 - (b[) 2 } = -4vHm(\;\ 2 7 ) . (23) 

Similarly, we can form another invariant from the vector b' and the vector appearing in the 
mass matrix: 

J 3 = kb'^ = M 12 b' 2 - M13&; = 2v / 2^ 3 Im(A 6 A;) . (24) 
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If we write the cubic interactions of Eq. [7| in terms of the physical scalars Xi in Eq. ^ as 
H-H+ ELi QiXi, we get 

J 2 = mlq k qi(T k2 T i3 - T k3 T i2 ) - mlT^qiqjT^T^T^ - T k3 T j2 ) , (25) 
J 3 = m 2 k T kl qi(T k2 T a - T k3 T i2 ) . (26) 

Thus, J 2 7^ or J 3 7^ imply CP violation. 

In the appendix we present a specific model in which there is no CP violation in the scalar 
mixing, but there is CP violation in the cubic and quartic scalar interactions. 

We are now in a position to discuss the cases in which S is a CP-odd monomial. We are 
looking for CP-odd invariants, and therefore, by carefully keeping track of the origin of the 
tensors we construct, we will know their CP properties and can use them to define invariants 
which are indeed CP odd. As an example we consider the cubic terms, 

2vG°H° [(lm\ 6 )R + (ReA 6 )/] , (27) 

This provides us with a CP-odd vector, hi = 2v (ImA6, ReAe) T = —eijbj/v, which we may com- 
bine with CP-even vectors or tensors such as fej and Cij to yield CP-odd invariants. Sometimes 
these invariants may be identically zero. For instance, 

- vbibi = Eijbibj = , (28) 
-vbiCijbj = biCije jk b k = J x . (29) 



3 More than two doublets 

We now extend the method of the previous section to the case in which there are more than 
two scalar doublets in the theory. Once again, we must work in the basis in which only one 
of the doublets, Hi, has a VEV v, and that VEV is real. Besides Hi there are now n more 
doublets in the theory, which may be written 

H " = ( (Ra +\)/V2 ) ' (30) 

with the index a going from 1 to n. We are free to choose a different basis for these n scalar 
doublets. That freedom is expressed by the invariance of the physics under the transformation 

H a — > K a pHp , (31) 

K being an arbitrary n x n unitary matrix. Let us write K as the sum of two n x n real 
matrices, Kr and Ki, in the following way: K = Kr + iKj. Let us moreover group all n fields 
R a and all n fields l a in a (2n)-dimensional vector u = (Ri, R 2 , R n , Ii, I 2 , I n ) T . Then, 
the effect of the change of basis in Eq. [31] on the fields R a and I a (we will consider later the 
charged fields) can be written u — > K^ 2n ^u, where the (2n) x (2n) matrix 

^ - ( % ~£ ) < 32 > 
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is orthogonal, as a consequence of K being unitary. Indeed, Eq. 
of the group U(n) in the group 0(2n). 

Similarly, a general CP transformation is written 



R a + il a -> U af3 (Rp - il ) 



constitutes an embedding 



(33) 
(34) 



with U a n x n unitary matrix. Separating U too in its real and its imaginary parts, U 
Ur + iUi, Eq. can be written u — > U^ n 'u, where 



% Uj 

Ui -u R 



(35) 



is a (2n) x (2n) orthogonal matrix. Notice however that is of a different form from K ^ . 

Let us now generalize the tensor of the previous section to a tensor e. This is a (2n) x (2n) 
matrix which is written in block form as 



e 








1„ being the n x n unit matrix. This matrix has the properties 

T 



U (2n) e \u (2n) 



e . 



-e . 



(36) 

(37) 
(38) 



The generalization of the results of the previous section is now evident. The coefficients of 
terms linear in the fields R a and I a are grouped in vectors, the coefficients of quadratic terms 
are grouped in second-rank tensors, and so on. Those quantities behave as vectors, second-rank 
tensors, and so on, under transformations of the sub-group U(n) of 0(2n) realized in Eq. |32[ 
When one wants to form a CP-odd quantity which is invariant under a basis transformation, 
one just picks some of these tensors and contracts all of their indices, leaving only two free 
indices % and j, which then are contracted with e^. 

As an example, consider, in the three-Higgs-doublet model, the mass terms of the neutral 
scalars, written in the form 



i ( H° R x R 2 h I 2 



M 



( H° \ 

Ri 

R2 
h 
h 



(39) 



M is now a 5 x 5 symmetric matrix. Besides the invariant M n , this matrix includes a vector 
hi = M l i+1 and a second-rank symmetric tensor = M i+1 j +1 . The tensor e can be written as 
a 4 x 4 antisymmetric matrix, with e 13 = e 24 = —£31 = —£42 = 1, and all other matrix elements 
0. The invariant J\ of the previous section can be extended to the three-Higgs-doublet model: 



J 1 — biCijbkB 



~-jk 



(40) 
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We now are able, however, to define other CP-odd quantities, for instance, 



J[ = hcijCjibkEik . (41) 

The analogous quantity for the two-doublet model (with e substituted by e) can be shown, with 
the help of Eq. to be equal to c&Ji, in which ca is itself a CP-even invariant. For more 
than two Higgs doublets, however, J[ is a CP-odd invariant independent of J\. In general, it is 
possible to construct still other independent CP-odd invariants, of ever increasing order in the 
elements of M. 

It is interesting to write the invariant J\ as a function of the masses of the five physical 
neutral scalars, and of the mixing matrix T, which is now 5x5. We obtain an obvious 
generalization of Eq. [19| 

Ji = m^rnjTiiTj^eabTi^+iTjt+i) . (42) 

It is important to note that it is the combination of T matrix elements in the right-hand-side 
of this equation that appears in the vertex of the Z° with the two physical particles X\ and Xj. 
Similarly, for J[ we find, 

J[ = m l 6 m|T a r i i(£ o6 Ti ia+ iT ii6+1 ) - (m 2 k Tl l )J 1 . (43) 



Eqs. and 43 actually hold for any number of Higgs doublets. 



4 Extension to the charged scalars 

The extension to the charged scalars of the work of the two previous sections is fairly obvious. 
As is seen in Eq. BTl the charged scalars H£ transform under basis transformations as the 



R a + il a . Together with Eq. g3|, this indicates that the procedure to follow must be to separate 
H£ into its real and imaginary parts, and to treat the real part as R a and the imaginary part 
as ./q, . 

For instance, the mass matrix of the charged scalars, N, is a n x n hermitian matrix, 
appearing in the mass terms 



( H ^ \ 



Hz ) N 



m 



(44) 



After separating each field H£ into its real and its imaginary part, and putting all the real 
parts, followed by all the imaginary parts, in a (2n) vector, these mass terms may be written 
using instead a (2n) x (2n) symmetric matrix which is, in block form, 



ReiV -ImiV 
imN ReiV 



(45) 



This matrix effectively constitutes a second-rank tensor which we may call hij. From this 
tensor, together with the tensor and the vector occurring in the mass matrix of the neutral 
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scalars, it is possible to construct further CP-odd quantities, like hijCj^b^bien, Cijhjkbkbi^ih or 
CijC jk h k ic lm e imj and so on. 

We emphasize that the treatment of the charged sector delineated in this section is necessary 
even in the two-Higgs-doublet model. Since the mass term for only one charged boson is of the 
form H~H + , and this is invariant under a rephasing of H + , the mass term of the charged scalar 
does not have any influence on CP violation. However, to deal with terms involving G~H + H° 
and G~H + £ il we must use our methods for the charged sector. The procedure consists in 
breaking into their real and imaginary component fields G R and Gj (G ± = G R ± iGj), and 
noting that they have CP even and CP odd properties, respectively. 

For example, in the cubic potential we have terms such as 



V2vG-H- 



+ h.c. 



Using = KeH + and £f = lmH + , we can rewrite this as, 

G R H°g& + G R ith& + • • • , 
where we have suppressed the analogous CP odd terms in Gj and used, 

V2, 



9i 



h 



V 

2V2 



v ■ - 2 

Therefore, we can recover J\ by considering transformations of the charged fields. 



(46) 

(47) 

(48) 
(49) 



5 Discussion and conclusions 

In this paper we addressed the problem of the systematic construction of CP-odd quantities 
for a SU(2)(g>U(l) model with many scalar doublets. We pointed out that basis trasformations 
leaving the VEV v invariant have no physical meaning. Therefore, our strategy has been to 
define CP-odd quantities which are invariant under those basis transformations. 

Our method does not discriminate between the two possible origins of CP violation, explicit 
or spontaneous. 

Our method does not rely on the assumption that the scalar potential is quartic. It can 
therefore be directly applied in effective scalar theories, with non-renormalizable interactions. 
The way is now open to search for experimental signatures of CP violation in the scalar 



sector |T3|, [10], and to determine which basic CP- violating quantities we might be able to 
measure. 



Appendix 

In this appendix we present a model with a symmetry such that there is no CP violation in the 
scalar mass terms, all the CP violation being in the cubic and quartic interactions. 
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Take three scalar doublets Hi, H 2 and H 3 , with a symmetry A (27) [|TT|, 0. This is the 
discrete subgroup of SU(3) generated by the two transformations 



Hi — > H 2 



H* 



Hi 



(50) 



and 



H x - Hi, H 2 - uH 2 , H 3 - cu 2 ^ , (51) 
where w = exp(27ri/3). The most general quartic Higgs potential invariant under this symmetry 



is 



V = fi(HlHi + H\H 2 + HlH 3 



+a 
+b 
+c 



\H\Hif + {H\H 2 f + (fljff s ) a j 
\h\Hi){hIh 2 ) + (HlHi)(HlH 3 ) + {HlH 2 ){H 
, „ L (HlH 2 )(HlHi) + {H\H 3 ){HtHi) + {H\H 3 ){HtH 2 ) 
+ {rf [{H\H 2 ){H\H 3 ) + {HlHi){H\H 3 ) + {HlHi){H\H 



+ /i.e.} 



(52) 



One possible vacuum of this potential has i?i with VEV t> while both i?2 and f/3 have 
vanishing VEV. Then, the only mass terms with a complex phase orig inate in d{H\H 2 ){H\H 3 ), 
because ci is complex. But it is obvious that we may redefine the phases of H 2 and H 3 in order 
to absorb the complexity of d in that mass term. On the other hand, that rephasing of H 2 and 
H 3 does not lead to the elimination of the complex phases in the terms d{H 2 Hi){H 2 H 3 ) and 
d(H 3 Hi)(H 3 H 2 ), which yield cubic and quartic interactions when Hi gets its VEV. Therefore, 
in this natural model there is no CP violation in the mixing of the scalars, but there is CP 
violation in their cubic and quartic interactions. 

Notice that this is what happens too in the model with A (27) symmetry of Branco, Grimus, 
and Gerard (BGG) ||11|| . There, the coupling d is real and CP is spontaneously broken, when 
Hi, H 2 and H 3 acquire VEVs with relative phases. However, in that model, when one goes to 
the basis in which only Hi has a non-vanishing VEV, one obtains a potential which still has 
A(27) symmetry, but now with a complex d. Therefore, the model of BGG has CP violation 
only in the cubic and quartic scalar interactions, but not in the mixing of the scalars. 
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